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ABSTRACT 


Recent data from the James Webb Space Telescope allow a first glimpse of galaxies at z = 11. 


^2 


The most successful tool for identifying ultra-high-redshift candidates and inferring their properties is 
photometric template fitting. However, current methods rely on templates derived from much lower- 
redshift conditions, including stellar populations older than the age of the Universe at z > 12, a stellar 
initial mass function which is physically disallowed at z > 6, and weaker emission lines than currently 
observed at z > 7.5. Here, two sets of synthetic templates, optimized for the expected astrophysics of 
galaxies at 8 « z « 12 and z > 12, are developed and used to fit three galaxies at z > 12 from the 
SMACS0723 field. Using these improved templates, quantitative estimates are produced of the bias in 
inferred properties from JWST observations at z > 8 due to these effects. The best-fit redshifts are 
similar to those found with previous template sets, but the inferred stellar masses drop by as much 
as 1-1.6 dex, so that stellar masses are no longer seemingly inconsistent with ACDM. The two new 
template sets are released in formats compatible with EAZY and LEPHARE. 


1. INTRODUCTION 


The launch of the James Webb Space Telescope 
(JWST) provides the first observational window into 
galaxies at z Z 11. Within a few weeks of the first 
data release, several galaxies have been found at z > 12 
(Adams et al. 2022; Atek et al. 2022; Castellano et al. 
2022; Donnan et al. 2022; Finkelstein et al. 2022; Naidu 
et al. 2022; Yan et al. 2022). The discovery of these 
galaxies and an analysis of their properties will allow 
tests of models for halo assembly and early galaxy evo- 
lution, and will likely be one of the key early results from 
JWST. 

Photometric surveys have been the best discovery tool 
for high-redshift galaxies (Steidel & Hamilton 1992), and 
all of the early, ultra-high-redshift sources have been 
detected in NIRCam imaging. Both the identification 
of these galaxies and the subsequent determination of 
properties such as stellar mass (M.) and star formation 
rates (SFR) rely on photometric template fitting. One 
of several codes (Brammer et al. 2008; Arnouts et al. 
1999; Leja et al. 2017; Carnall et al. 2018) is used to 
compare a series of model spectral templates against 
the observed photometry, and the redshift and proper- 
ties are inferred from the best-fit models. Most models 
are synthetic templates produced from stellar popula- 
tion synthesis under a range of assumptions and con- 
ditions. Crucially, all these synthetic templates are de- 
rived from local conditions, which are either theoreti- 


cally disallowed or observationally disproven for galaxies 
at ultra-high redshift. 

Therefore, in this work, a new template set is de- 
veloped with the aim of producing models tuned for 
fitting ultra-high-redshift galaxies. These templates 
are also available for download at https://github.com/ 
e-m-garcia/hot-templates in formats compatible with 
EAZY and LEPHARE. The models are significantly dif- 
ferent in three respects. Firstly, they use younger stellar 
populations since at those redshifts, the stellar popula- 
tion cannot be more than a few hundred Myr old. In 
contrast, all current templates are older than the age 
of the universe at z ~ 12 (Brammer et al. 2008). Sec- 
ondly, we apply bottom-lighter IMFs corresponding to 
expected conditions in star-forming regions of galaxies 
at z ~ 15 (see 8 2). Thirdly, we use stronger emission 
lines as none of the local strong-line metallicity calibra- 
tions provide a good prediction of the observed metal- 
licities at z > 7.5 (Curti et al. 2022). The differences in 
inferred parameters between these templates and ones 
derived from current templates are discussed in 8 3. 

Throughout this paper we assume a flat ACDM cos- 
mology with €,,9 = 0.3, Qa,o = 0.7 and Ho = 70 km 
s-1 Mpc-!. 

2. THE HIGH-REDSHIFT INITIAL MASS 
FUNCTION 


The physical parameters inferred for high-redshift 
galaxies will depend upon the assumptions made when 
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producing the synthetic spectra used to fit them. Under 
the correct assumptions, templates should produce both 
a better fit and a more accurate description of the key 
physical properties of the galaxy than with incorrect as- 
sumptions. One central assumption common to all cur- 
rent templates is a universal stellar initial mass function 
(IMF), which is assumed to be identical to the current 
Galactic IMF. This is a critical assumption for infer- 
ring properties, because the light emitted by a galaxy is 
dominated by the most massive stars comprising only a 
small fraction of the stellar mass, and the IMF is nec- 
essary to infer the remainder of the stellar population 
from that high-mass tail. However, theoretical mod- 
els for star formation predict that the IMF should de- 
pend upon several properties of star-forming molecular 
clouds, including gas temperature and metallicity (Low 
& Lynden-Bell 1976; Larson 1985; Jermyn et al. 2018; 
Steinhardt et al. 2020). Thus, assuming a universal IMF 
is, in effect, an assumption that the differences between 
Galactic and extragalactic star-forming regions are neg- 
ligible. 

The ideal approach would be to construct templates 
using a range of possible initial mass functions and 
determine which one fits best. This is the approach 
taken in recent work fitting galaxies at lower redshifts 
(Sneppen et al. 2022; Steinhardt et al. 2022a,b). Based 
on theoretical arguments (Jermyn et al. 2018), a one- 
parameter family of templates at different gas tempera- 
tures 7; was produced. The corresponding IMF at each 
gas temperature was taken as 


2 
m 93 m <0.08Mo - (&) 


which matches a Galactic (Kroupa) IMF at T, = Tp = 
20 K. Other investigations of the relation between tem- 
perature and IMF produce a similar relationship, but 
with a slightly different temperature dependence (Jeans 
1902; Hopkins 2012; Chabrier et al. 2014). Each of 
these sets of templates was used to fit photometric ob- 
servations from the COSMOS2015 catalog (Laigle et al. 
2016), and the best-fit among all of the template sets 
was taken as having the correct IMF and used to infer 
physical parameters. 

However, as there are strong spectral covariance be- 
tween changes in the IMF and changes in dust or metal- 
licity, it was only possible to constrain the IMF for the 
~ 1096 of the catalog with the highest signal-to-noise 


ratio in Sneppen et al. (2022). Further, the COSMOS 
catalog includes multi-wavelength measurements in up 
to 30 bands. Thus, a similar approach used on NIR- 
Cam photometry, with typically 4-5 bands exhibiting 
non-zero flux emitted by the most distant observable 
galaxies, it is not possible to observationally constrain 
the IMF. Rather, it will be necessary to simply choose 
an IMF (or, equivalently, extinction and metallicity) and 
assume it holds for all galaxies in the sample. Indeed, 
this is already being done by all template fitting codes, 
which assume that the high-redshift IMF is the same as 
the Galactic one. 

An alternative might be to produce observational 
rather than synthetic templates. Using NIRSpec, it 
should be possible to take spectra of a number of these 
z ~ 15 sources, and those with confirmed high redshifts 
can then be used as templates to fit others. A simi- 
lar empirical approach is taken for one of the template 
sets in LEPHARE (Arnouts et al. 1999). However, the 
spectrum alone will not be able to constrain the IMF, 
and thus this approach cannot be used to infer proper- 
ties such as stellar mass or star formation rate. If the 
sole goal is redshift determination, this may still be the 
best source of templates, as the templates will certainly 
correspond to at least some types of z ~ 15 galaxies. 

'The remaining approach, then, is to simply assume 
az ~ 15 IMF and hope that the assumption is ap- 
proximately correct. This is already being done with 
the assumption of a Galactic IMF, and the goal here 
is to improve upon that assumption. There are three 
key observational and theoretical indications validating 
a change of IMF in this ultra-high-redshift regime: 

First, the best-fit IMFs for COSMOS galaxies exhibit 
several clear trends (Sneppen et al. 2022). Quiescent 
galaxies at all redshifts out to z ^ 2 are well fit with 
approximately-Galactic IMFs (Steinhardt et al. 2022b). 
However, for star-forming galaxies, there is a charac- 
teristic, best-fit IMF at each redshift, which becomes 
bottom-lighter towards high redshift (Steinhardt et al. 
2022a). Inferred values of T, rise from 25 K at low 
redshift to over 30 K by z — 4. An extrapolation of 
this characteristic T; to z ~ 15 produces an estimate of 
35 — 40 K. 

Second, a reasonable interpretation is that in star- 
forming galaxies, the young stellar population is the 
dominant contribution to gas temperature in star- 
forming regions. If so, z ~ 15 galaxies, which must have 
very high SFR in order to have formed so quickly, should 
have even higher T} and even bottom-lighter IMFs. A 
closer examination of the first stages of star formation, 
taking place primarily in galactic centers, finds T, of 


45 — 60 K (Steinhardt et al., in prep.), with 60 K char- 
acteristic of galaxies in their earliest stages of evolution. 

Finally, at high redshift the CMB temperature is an 
additional source of heating. At z ~ 15, the CMB tem- 
perature is 44 K, and if those stars have been forming 
over the past ~ 120Myr, a typical age for galaxies on 
the star-forming main sequence, a luminosity-weighted 
CMB temperature at the time of star formation would 
be closer to 55 K. It is difficult to find mechanisms for 
star-forming clouds to cool well below the CMB temper- 
ature, so this should be a minimum possible value of T; 
at high redshift. 

All three of these estimates suggest that a reasonable 
value of T, for galaxies at z ~ 10 would be 40 — 50 K 
and at z ~ 15 significantly higher. In this work, the 
choice is made to develop templates based on T, = 60 
K for use when modeling the highest-redshift galaxies in 
JWST at z > 12. This produces a significantly bottom- 
lighter IMF than in the Milky Way. An additional set of 
templates with T; = 45 K is also provided, and is likely 
a better model for galaxies at 8 < z « 12. 

In summary, inferring properties from JWST observa- 
tions creates somewhat of a quandary. Ultra-high red- 
shift galaxies with very young stellar population are si- 
multaneously (1) the galaxies for which inferred prop- 
erties are most sensitive to the choice of IMF; (2) the 
galaxies most likely to form stars under very different 
conditions than in our own Galaxy today; and (3) the 
galaxies for which it will be most difficult, if not impos- 
sible, to constrain the IMF observationally. 

Fortunately, measurements from large surveys at 
lower redshift do provide constraints on the IMF, and 
their measured IMFs are consistent with physical intu- 
ition. Thus, the choice made in this work is to use sim- 
ilar physical intuition to extrapolate to the most likely 
ultra-high-redshift IMF, and this is almost certainly a 
better approximation to the true IMF than a Galactic 
one. However, it would not be shocking to discover that 
this physical intuition breaks down under the extreme 
conditions in the earliest galaxies, and inferred proper- 
ties such as stellar masses and star formation rates for 
ultra-high-redshift galaxies should be assumed to have 
systematic uncertainties well in excess of those produced 
using the same techniques at lower redshift. 


3. FITS AND INFERRED PARAMETERS 


'The IMFs described in this work are implemented for 
two different photometric template fitting codes, EAZY 
(Brammer et al. 2008) and LEPHARE (Arnouts et al. 
1999). EAZY starts from a set of 560 model spectra 
generated using Flexible Stellar Population Synthesis 
(FSPS), which span a wide range of SFR, stellar popula- 
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tion ages, metallicities, and extinctions, then condenses 
them into 12 basis vectors which approximately span 
the full space!. Within that limited basis, matrix inver- 
sion allows a rapid determination of the best-fit linear 
combination of those 12 templates to the observed pho- 
tometry over a grid of redshifts. LEPHARE, instead, 
relies on a grid search. For a smaller number of star 
formation histories, a grid of ages, redshifts, and extinc- 
tions is compared against the observed photometry to 
produce the best fit. LEPHARE will therefore only con- 
sider physically reasonable models, although the grid 
search means that often only a coarse search is possible 
in a reasonable runtime and the best-fit model might 
not be found. EAZY can find the best fit more quickly 
within the limited basis, but in principle the best-fit lin- 
ear combination might be physically unreasonable if the 
true parameters lie well outside of that basis. 
Implementation for a grid-based search such as LEP- 
HARE is straightforward, as it merely requires generating 
templates for every combination under consideration at 
a variety of ages. Instructions for installing the template 
set in this work for use with EAZY (or other codes that 
use the same format, e.g., STARDUST; Kokorev et al. 
2021) and LEPHARE can be found in Appendices A and 
B, respectively. The remainder of this section focuses 
on the template sets produced for EAZY, which addi- 
tionally require a reconsideration of the template basis. 


3.1. Modifying the Standard Template Basis 


The standard EAZY release includes a set of 12 basis 
templates with properties chosen to approximately span 
the full observed galaxy population in large photometric 
catalogs. A list of the templates and their properties can 
be found in Table 1. The standard EAZY templates are 
based on a Chabrier IMF, but the IMFs used here (Eq. 
1) are instead equivalent to a Kroupa IMF at T, = 20 K. 
Therefore, a set of templates using a standard Kroupa 
IMF is used for comparison in order to demonstrate the 
effects of varying the IMF on inferred properties. 

Each template is a specific linear combination of a 
broader set of 560 individual models. To produce tem- 
plates at T = 45 K, each of the 560 individual models 
was run with a modified IMF as in Eq. (1). Those 
templates were then combined into the 12 basis tem- 
plates using the same coefficients. This produces a very 
similar set of luminosity-weighted average properties for 
each template (see Table 2). Because of the bottom- 
lighter IMF, the light in each template is dominated by 


1 For some catalogs, additional templates are added to produce an 


expanded basis. 
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Table 1. Properties of the stellar populations in the 12- 
template basis of EAZY. The properties are nearly identical 
for different IMFs, as they are constructed from linear combi- 
nations of FSPS models with the same physical parameters. 
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Age (Gyr) | Av Z 
0.72 0.06 | 0.0066 
3.15 0.11 | 0.0119 
10.59 0.13 | 0.0140 
13.83 0.12 | 0.0223 
13.64 0.73 | 0.0237 
15.12 2.88 | 0.0240 
0.34 0.07 | 0.0154 
0.62 0.27 | 0.0211 
1.07 0.80 | 0.0216 
1.30 1.37 | 0.0216 

10 2.26 1.81 | 0.0217 

11 2.44 2.94 | 0.0216 
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slightly younger stars which have produced a slightly 
higher extinction and metallicity than at T = 20 K. 

For illustration, these templates have been used to fit 
three high-z candidates (z > 12) from the SMACS0723 
cluster field (Kokorev et al. in prep.). These have been 
labelled as SMACS 1, 2 and 3. The data reduction of 
the raw images and the photometry extraction have been 
performed with the GRIZLI pipeline (Brammer & Math- 
aru 2021). The images have been additionally aligned 
to match the astrometric reference frame used for the 
ALMA Lensing Cluster Survey data reduction of HST 
data (Kokorev et al. 2022). The data reduction and the 
photometric catalog will be fully described in Brammer 
et al. 2022 (in prep.). The best-fit redshifts using the 
standard EAZY templates are z = 12.23, 13.48 and 
15.09, respectively. The 60 K Kroupa IMF templates 
produce a very similar redshift, which is not surpris- 
ing given that these are dropout galaxies in the JWST 
FO90W band and the redshift is primarily determined 
from the presumed location of the break. The fits them- 
selves are also broadly similar, with typical differences 
of < 2% in the reconstructed flux (Fig. 1). 

The physical properties exhibit a far greater difference 
(see Fig. 2 or Table 2). In a star-forming galaxy, most 
of the light comes from the most massive stars, while 
most of the stellar mass lies in the least massive ones. 
Thus, a significantly bottom-lighter IMF will produce a 
much lower stellar mass and SFR. 

The best-fit stellar masses are between 0.5 and 0.6 dex 
lower at T, = 60 K than at T, = 20 K. At lower red- 
shifts, differences in the best-fit extinction and metallic- 
ity create a more complex relationship between T, and 
the best-fit stellar mass (Sneppen et al. 2022; Steinhardt 
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Figure 1. Best-fit reconstructed spectra (left) and red- 
shift probability functions (right) for three galaxies in 
SMACS0723. The fits were performed with the standard 
EAZY Chabrier IMF at T, = 20K (blue) and with Kroupa 
IMF at T; = 60K (maroon), with the latter fit expected to 
be a better description of galaxies at z > 12. 


et al. 2022b). However, galaxies at z ~ 15 have not had 
enough time to form as much dust. Thus, this change 
in stellar mass is very close to the difference in mass- 
to-light ratios for a single stellar population between 
T, = 20 K and 60 K (Fig. 2). A similar shift exists 
in the best-fit star formation rates. 


3.2. Younger Stellar Templates for High-Redshift 
Galaxies 


With either set of fits, the inferred specific star for- 
mation rates (SSFR, or SFR/M,) are ~ 10 5/yr?. At 
that rate, the entire stellar population would form in 
100 Myr. At z — 12.4, 14.3, and 15.7, the Universe is 
350 Myr, 290 Myr, and 260 Myr old, respectively. Thus, 
100 Myr is a plausible age for a stellar population in a 
z~ 15 galaxy. 

However, every one of the basis templates corresponds 
to an older stellar population than in these high-redshift 


? In the standard EAZY basis, a higher sSFR cannot be inferred, 
since none of the 12 basis templates has an sSFR over ~ 10-8 /yr. 


Table 2. Physical properties for the best-fit template reconstruction of three galaxies in the SMACS0723 
field. Four template sets are compared: (1) the standard EAZY basis, derived using a Chabrier IMF; (2) 
the same basis derived with a Kroupa IMF; (3) the same basis at T; = 60K; and (4) a new T; = 60 basis 
with younger stellar populations and other properties tuned for ultra-high-redshift galaxies. All template 
bases agree that these sources lie at z > 12, but the inferred masses differ far more significantly. There is 
insufficient information to use a goodness-of-fit test to determine which of these template sets will fit best, 
so the approach recommended here instead relies on astrophysical assumptions about these galaxies. 


IMF Template Set log,j((M/Mc) SFR (Mo yr ) t (Gyr) Av z 

SMACS 1 (07:22:44.85, -73:29:53.75) 

Chabrier (20 K) EAZY 12 8.851000 ° 3.02 0.34 0.07 12.23 

Kroupa (20 K) EAZY 12 8.64* 0-05 4.28 0.38 0.07 12.41 

Kroupa (60K) ^ EAZY 12 813 02 2.45 0.39 0.07 1241 

Kroupa (60 K) ^ Modified 6 BET 2.13 0.00 0.10 12.06 
SMACS 2 (07:22:51.93, -73:29:21.89) 

Chabrier (20 K) ^ EAZY 12 8.911000 3.47 0.34 0.07 13.48 

Kroupa (20K) | EAZY 12 8.604062 3.32 0.34 0.07 14.26 

Kroupa (60K) ^ EAZY 12 8.064062 1.54 0.33 0.07 14.30 

Kroupa (60 K) Modified 6 Bde 1.78 0.13 0.13 13.69 
SMACS 3 (07:23:03.45, -73:28:47.01) 

Chabrier (20 K) ^ EAZY 12 9.00* 0-06 4.29 0.34 0.07 15.09 

Kroupa (20K) | EAZY 12 8T 5 cs 4.25 0.34 0.07 15.77 

Kroupa (60K) ^ EAZY 12 8.167005 1.97 0.33 0.07 15.80 

Kroupa (60 K) ^ Modified 6 caer 1.72 0.10 0.01 15.23 


? Uncertainties shown here are those reported from the EAZY covariance estimation. As described in 
§ 3.2, these fit uncertainties significantly underestimate the true uncertainties for reconstructions consist- 
ing only of one template. Thus, negligible uncertainties do not indicate a very well-constrained fit, but 
rather that the true SED lies well outside of the template basis. 


galaxies. The template with the youngest stellar popu- 
lation age, template 6 (see Table 1), has a luminosity- 
weighted age of 340 Myr. Nevertheless, as the youngest 
stellar population, it is the closest match. The best-fit 
linear combinations of the 12 EAZY templates are dom- 
inated by this single youngest template, with template 
6 contributing at least ~ 76% to the fits of the object 
1 and 100% to the other fits. The next-youngest tem- 
plate, template 7, contributes all of the remainder. The 
true SED very likely lies outside of the basis spanned by 
these 12 templates. 

If the true SED lies well outside of the basis spanned 
by these 12 templates, the closest match will likely lie 
at an extreme edge of that basis. This is consistent with 
the EAZY best-fit reconstructions, which consist almost 
solely of a single template. If measurement uncertainties 
perturb the input photometry, it will still result in a 
galaxy younger than any of the basis templates, and so 
the best fit will still be dominated by the same single 
template. Thus, even a large uncertainty in photometry 
will result in a minimal change in the best-fit spectrum 
and parameters. 


Because the uncertainties on EAZY fit parameters are 
derived from the uncertainty on the best fit within this 
12-template basis, the estimated uncertainty on param- 
eters such as stellar mass will also be unreasonably low. 
The EAZY uncertainty estimation in the stellar masses 
for the three SMACS galaxies using the standard tem- 
plate basis breaks, as the covariance matrix cannot be 
calculated. The resulting reported zero uncertainties are 
not an indication of reliability in the masses, but rather 
of how far outside of the template basis the observed 
photometry apparently lies. 

Here, a set of six additional templates is introduced to 
solve this problem. They include three different stellar 
population ages and two choices of extinction at each 
age, all of the same metallicity (Table 3). Two of these 
templates also have increased equivalent widths of emis- 
sion lines to EW(Hf--[O 11]\5007) ~ 1370 A. Such ex- 
treme nebular emission, thought to be associated with 
the “IRAC excess” in galaxies at z > 7, was observed 
in IRAC photometry (eg., Smit et al. 2014; Oesch et al. 
2015; Roberts-Borsani et al. 2016), NIRCam photome- 
try (Labbe et al. 2022), as well as NIRSpec spectra (Car- 
nall et al. 2022). The six extended templates T, = 45 K 
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Table 3. Properties of the stellar populations in the 6 new 
FSPS templates for EAZY at T, = 45 K and 60 K. The tem- 
plates all have young stellar populations with ages between 
0.05 and 0.20 Gyr for 45 K and between 0.1 and 0.5 Gyr for 
60 K. The low extinction and metallicity are expected to be 
typical of ultra-high-redshift galaxies. Because of the strong 
emission lines present in early NIRSpec observations (Car- 
nall et al. 2022), two templates are included with strongly 
enhanced EW(Hg + [O 11]\5007) with respect to the stan- 
dard FSPS setting of ~ 274 A. 


ID | tas (Gyr) teo (Gyr) Av,45,60 245,60 Neb. Em. 
1 0.1 0.05 0.005 0.0012 | 1x (std.) 
2 0.1 0.05 0.5 0.0012 5x 
3 0.3 0.1 0.005 0.0012 | 1x (std.) 
4 0.3 0.1 0.5 0.0012 | 1x (std.) 
5 0.5 0.2 0.005 0.0012 | 1x (std.) 
6 0.5 0.2 0.5 0.0012 5x 


and 60 K are available online, with installation instruc- 
tions in Appendix A. 

'The two different values for extinction of the new tem- 
plates are Ay — 0.005 and Ay — 0.5. Template 6, the 
dominant template in the fits with the standard EAZY 
basis with all three IMFs considered, has Ay = 0.06, 
which is therefore the extinction fit for all three galax- 
ies using those bases. Although the fits with the new 
six-template basis could produce extinction as high as 
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Figure 2. Expected change in inferred stellar masses using 
the IMFs in this work relative to Galactic IMF as a function 
of Tg. The dashed black line indicates the average shift based 
on the analytical integral while holding all other parameters 
fixed, which is a close match for the observed differences in 
fitting the three z > 12 candidates under the same condi- 
tions. The actual shift in stellar mass will combine this ef- 
fect, the difference between Kroupa and Chabrier choices of 
Galactic IMF, and changes in stellar age and nebular emis- 
sion lines and will typically exceed the analytical integral 
shown here. 


Ay — 0.5, all three SMACS objects are fit with much 
lower extinctions similar to Template 6. 

Although large multi-wavelength surveys can have as 
many as 30 photometric bands, which can be parameter- 
ized by a basis of 12 templates, these ultra-high-redshift 
candidates will typically will only have detections in 4- 
5 bands. Thus, the basis we advocate using for JWST 
objects is necessarily smaller than the full EAZY basis. 

Finally, the ability of the new templates to accurately 
constrain redshift was cross-checked by fitting an object 
in the relevant redshift range detected spectroscopically 
with NIRSpec, which was first reported in Carnall et al. 
(2022) (object ID 4590 in the paper). The best-fit red- 
shift at Ty = 45 K is Zphot = 8.32, which agrees with 
Zspec = 8.498 well within one ø of the zj5;; distribution. 


4. DISCUSSION 


'This paper presents a revised and improved basis of 
templates for fitting high-redshift galaxies with z > 8. 
While current synthetic templates are derived from stel- 
lar populations, which are impossibly old, have weaker 
emission lines than observations and invoke an unphysi- 
cal assumption of a Galactic IMF, the revised templates 
improve on all these accounts. Here, a new set of six 
templates optimized for the extreme conditions in the 
earliest galaxies is presented, along with instructions for 
use with EAZY and LEPHARE. 

It should be noted that the addition of these templates 
to the standard 12 EAZY basis templates will provide 
additional degrees of freedom, and JWST observations 
alone will typically not provide enough information to 
constrain the best fit within a combined, 18-template 
basis. Further, the inclusion of very young stellar popu- 
lations might allow spurious fits for some lower-redshift 
galaxies, since in combination with dusty or quiescent 
templates, it allows matching the blue and red ends of 
the SED nearly independently. For similar reasons, de- 
riving a best-fit IMF for each individual source is not 
possible, so it is instead necessary to select a single set 
of templates based on the expected astrophysics of typ- 
ical very high-redshift galaxies. 

Therefore, to avoid improving high-z fits and prop- 
erties at the cost of degrading low-z objects we recom- 
mend a two-stage procedure. First, the standard EAZY 
template basis should be used to fit all sources, since it 
is well optimized for the local objects in a photometric 
survey, which constituted the vast majority. Then, the 
objects with best-fit z > 8 should be refit with these new 
templates: the template set at T; = 45 K at 8 < z < 12 
and with Tj, = 60K at z > 12. The latter fits should 
produce better estimates of the physical parameters for 
these galaxies, which are robustly fitted at the highest 


redshifts. The new templates introduced are intended as 
a first-order improvement to current low-redshift derived 
templates, but are certainly not a complete solution. Fu- 
ture spectroscopy at z>12 may provide improved prop- 
erties for template-fitting at these redshifts. 


4.1. Effects on Redshift Determination 


For galaxies at z < 4 in COSMOS, a change in IMF 
had negligible effect on the best-fit redshifts (Sneppen 
et al. 2022). Most objects have detected flux on both 
sides of the Balmer break, often from narrow bands, 
and therefore the location of the break was well con- 
strained. However, the highest-redshift JWST galaxies 
are instead found as dropouts in NIRCam photometry, 
with no narrow bands to further constrain the location 
of the break. The bluest bands will have no detection, 
the reddest bands will have clear detections, and one 
band in between (F200W at z ~ 15) will typically have a 
fainter detection, combining part of the spectrum above 
the Lyman break with negligible flux at shorter wave- 
lengths. The strongest constraint on the best-fit redshift 
comes from essentially calculating the equivalent width 
for that partial detection and determining what fraction 
of this band has non-zero flux. Thus, a change in the UV 
slope near the Lyman break will change that equivalent 
width and unlike at z < 4, the corresponding best-fit 
redshift will be sensitive to the IMF. 

For example, with everything else held constant,the 
switch from a Chabrier IMF to a Kroupa IMF within 
the standard EAZY template basis yields a higher best- 
fit redshift. Because the best-fit models with this 12- 
template basis are solely composed of the same single 
template with both IMFs, a higher redshift is found for 
each SMACS source. These changes in redshift are still 
within the estimated redshift probability density func- 
tion, which describes the statistical uncertainty. How- 
ever, they are not statistical, but rather a systematic 
offset between the redshifts calculated from each tem- 
plate set due to smaller equivalent widths with a Kroupa 
IMF. 

With the six-template basis proposed in this work, a 
more complex dependence is possible. In addition to a 
change in IMF, the additional templates allow galaxies 
to be fit with different metallicities and extinction values 
than using the standard EAZY basis. Thus, in total the 
UV slope of the reconstructed spectrum could be either 
steeper or shallower, and thus the redshift could change 
in either direction. 


4.2. Impossibly Early Galaxies 


The strongest effect will be on the inferred stellar 
masses. The change between a Chabrier and Kroupa 
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IMF can produce a difference of 0.2 — 0.3 dex in stellar 
mass, but the effect of a high Ty is far stronger with 
variations of ~ 0.5 — 1 dex. 

This change is particularly relevant at very high red- 
shift, because the inferred stellar masses from standard 
template libraries are high enough to challenge the stan- 
dard ACDM paradigm. Prior to JWST, the most mas- 
sive, highest-redshift galaxies in large photometric sur- 
veys were already difficult to reconcile with ACDM halo 
mass functions (Steinhardt et al. 2016; Behroozi & Silk 
2018), an effect described as the impossibly early galaxy 
problem (Steinhardt et al. 2016). One possible solu- 
tion would be an increase in the stellar baryon fraction 
(Finkelstein et al. 2015; Behroozi et al. 2020), which 
would allow halos of the same mass to contain galax- 
ies with greater stellar masses at high redshift. How- 
ever, the estimated masses of JWST galaxies at z = 10 
(Labbe et al. 2022) are high enough that even this is 
no longer a solution: if the inferred stellar masses and 
redshifts are correct, the stellar masses would exceed 
the full baryon masses within halos at z = 10 (Boylan- 
Kolchin 2022). 

The T, = 45 K templates should provide better fits 
for z = 10 galaxies than templates with a Galactic IMF. 
However, these templates cannot satisfactorily match 
the observed photometry, since these z = 10 candidates 
were selected because they appear to exhibit strong 
Balmer breaks. None of the templates in the 6-template 
basis have sufficiently strong Balmer breaks, which take 
Z, 300 Myr to develop for a single stellar population 
and = 500 Myr for any of the more realistic star forma- 
tion histories used in these templates. The age of the 
Universe at z — 10 is less than 500 Myr. Labbe et al. 
(2022) solve this problem by using an older and dustier 
stellar population than the T; = 45 K template set to 
produce a strong Balmer break. This reconstructed stel- 
lar population is significantly older than 500 Myr and 
thus physically impossible. Given current models, it is 
not possible to produce as strong of a Balmer break at 
z = 10 as exists in the reported photometry. 

Perhaps the most likely explanation is that im- 
proved zero point corrections will significantly reduce 
the strength of the break. Another possibility is that 
these are very dusty, much lower-redshift sources. Con- 
straining the redshifts to assume these are truly z ~ 10 
sources and using the 7; = 45 K templates, the stel- 
lar masses for typical z — 10 galaxies drop by 1.6 dex. 
'This large difference comes from a combination of sev- 
eral effects: (1) a bottom-ligher IMF; (2) the difference 
between Kroupa and Chabrier IMFs; (3) a younger stel- 
lar population; (4) lower extinction; and (5) enhanced 
nebular emission lines producing a fainter continuum. 


This last effect could provide an alternative explana- 
tion for the apparent Balmer breaks in these z = 10 
galaxies. NIRSpec observations find stronger emission 
lines at z > 7.5 (Carnall et al. 2022; Harikane et al. 2022) 
than are present in FSPS templates (Conroy et al. 2009; 
Conroy & Gunn 2010) derived from lower-redshift ob- 
servations. Sufficiently strong emission lines might even 
mimic a spectral break. If so, followup spectroscopy of 
these massive z ~ 10 candidates might instead find very 
strong emission lines but no discernible Balmer break. 
This would be consistent both with theoretical expecta- 
tions for a stellar population when the Universe is less 
than 500 Myr old and with the observed combination 
of strong emission lines without Balmer breaks in the 
Carnall et al. (2022) NIRSpec spectroscopy. 

Although the Galactic IMF-derived stellar masses are 
inconsistent with ACDM, the T, = 45 K masses can 
still be reconciled with theoretical mass functions (Fig. 
3). More generally, the high-luminosity tail of the lu- 
minosity function at very high redshift would no longer 
necessarily correspond to stellar masses as high as at 
more moderate redshifts. As selection improves, it is 
possible that additional massive objects will be discov- 
ered, and the stellar mass densities estimated at present 


might be underestimated. If so, a complete survey might 
still find tension between theoretical models for halo for- 
mation and observed high-redshift galaxies. However, 
at present the impossibly early galaxy problem can be 
solved solely with a change to a template basis consis- 
tent with CMB heating. 

Spectroscopy and improved modeling will ultimately 
be needed to determine whether these galaxies are truly 
at the stellar masses inferred from template fitting with 
any of these bases. Currently, however, one must as- 
sume an IMF in order to determine physical parame- 
ters, which are essential in order to use JWST observa- 
tions to test fundamental astrophysics and cosmology. 
Given the extreme conditions at very high redshifts, it 
is necessary to use templates derived from a physically- 
motivated IMF rather than one derived from local obser- 
vations under conditions that are physically impossible 
at high redshift. 
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A. INSTALLING TEMPLATES IN EAZY 


EAZY templates constructed for high redshift studies consist of two two sets of 6 templates. Each set has the IMF 
corresponding to gas temperatures T} = 45 and 60 K. They were produced using FSPS (Conroy et al. 2009; Conroy & 
Gunn 2010), with the 6 new templates representing stellar populations younger, more metal poor than in the standard 
12 basis EAZY set and with stronger nebular emission (see Table 3). 

For more detailed information on usage and further resources, refer to the Github repository at https://github.com/ 
e-m-garcia/hot-templates. EAZY with the standard set of templates is available at https://github.com/gbrammer/ 


eazy-py. 


B. INSTALLING TEMPLATES IN LEPHARE 


'The set of templates produced for usage consist of 4 templates at 45 K and 60 K each of varying star-formation 
histories at a z—0.0012 metallicity. These are reconstructed in FSPS (Conroy et al. 2009; Conroy & Gunn 2010) which 
allows for replacing the original Chabrier IMF (Chabrier 2003) with a temperature-dependent Kroupa IMF (Kroupa 
2001; Jermyn et al. 2018). The templates are written into a BC03 ASCII file format such that LEPHARE can read 
them. 

'To use the templates, simply include the given templates in the galaxy template directory in LEPHARE, and use the 
young stellar age list included. It is also recommended to start from parameters associated with the existing 12 BC03 
templates used for COSMOS when reading the templates, constructing the magnitude grid, and estimating redshifts 
and properties, although these can be modified to fit the needs of the fitting (such as different filters). For more detailed 
information on usage and further resources, refer to the Github at https://github.com/e-m-garcia/hot-templates and 
the LEPHARE repository at https://gitlab.lam.fr/Galaxies/LEPHARE. 


